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Abstract 

We calculate the van der Waals dispersive interaction between a neutral but polarizable atom and a 
perfectly conducting isolated sphere in the nonretarded case. We make use of two separate models, one 
being the semiclassical fluctuating-dipoles method, the other using ordinary quantum mechanics. 



1 Introduction 

The existence of attractive intermolecular forces be- 
tween molecules of whatever kind has been long as- 
sumed due to the possibility of every gas to liquefy. 
The first quantitative (although indirect) characteri- 
zation of these forces was done by J.D. van der Waals 
in his 1873 thesis [1] in the equation of state for real 
gases, which can be written as 



P 



—) {V-b) = nRT 



(1) 



in which P, V, T and n are the pressure, volume, 
temperature, and the number of mols of the gas un- 
der consideration, respectively, and R is the universal 
gas constant. The parameters a and b - the van der 
Waals constants, which vary from one gas to another 
- can be evaluated through fitting this equation with 
experimental data. While the parameter b relates to 
the finite molecule volume, and to the fact that an 
individual molecule cannot access the entire volume 
V of the gas, the term {a/V"^) relates to an attrac- 
tive intermolecular force. Such forces received the 
general name of van der Waals forces. It should be 
stated that when one mentions intermolecular forces 
as such, one assumes that the separation between 
the molecules in question (or atoms, in the case of 
monoatomic gases) is large enough as to exclude the 
overlapping of electronic orbitals. These forces are 
usually distinguished in three types, to wit: orienta- 



tion, induction and dispersion van der Waals forces, 
which we briefly discuss - more detailed discussions 
can be found in references [2], [2], [4]. 

Orientation forces occur between two polar 
molecules, i.e. two molecules possessing per- 
manent electric dipoles, e.g. water molecules. 
These forces were first computed by W.H.Keesom, 
in the early twenties, considering the ther- 
mal average of the interaction energy of two 
randomly-oriented electric dipoles di and d2, 
namely, [-2/?,){d\dl/ATTeor^){l/kBT) for keT > 
di'^2/47reor^, where r is the distance between the 
dipoles (molecules) and fcs is the Boltzmann con- 
stant. Although the amount of possible attractive 
orientations equals the amount of repulsive ones, once 
we take into account that attraction setups corre- 
spond to smaller energies and that the Boltzmann 
weight is ^-^I^bT jg^ diminishes with in- 

creasing energy £"), it can be easily understood why 
orientation forces are attractive. We also note that 
orientation forces decrease with increasing tempera- 
ture, which is natural since higher temperatures turn 
repulsive orientations as accessible as attractive ones. 

It was recognized by P.Debye and others that 
there ought to exist an interaction between a polar 
molecule and an apolar, albeit polarizable, one, once 
the polar molecule induces a dipole in the other one, 
giving rise to a dipole-dipole attraction force, respon- 
sible for the induction van der Waals forces. (In fact. 
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even a quadrupole or higher permanent multipole can 
induce such a dipole and give rise to induction forces.) 
The attractive character of these forces foUow from 
the fact that an induced dipole is parallel to the in- 
ducing field (in the case of non-isotropically polar- 
izable molecules this is at least approximately true) 
and that the interaction energy between an electric 
field and a parallel (or close to parallel) dipole is al- 
ways negative. This correlation leads to a nonvan- 
ishing force at increasing temperatures. Moreover, 
we can evaluate the dependence of the force on the 
distance r between the molecules recalling that the 
magnitude of the field Ei generated by the perma- 
nent dipole di is proportional to di/r^ , and that the 
energy of the interaction between this field and the 
second (induced) dipole (d2) is of the form 



U{r) 



-aEl{r) oc 



(2) 



where a is the molecular polarizability (assumed lin- 
ear, for simplicity). Since U{r) is proportional to 
(— 1/r^), the force has a (— 1/r'') behavior. 

It just so happens that the correction term [a/V"^) 
can enhance the ideal gas approximation for each and 
every gas known in nature, including gases consti- 
tuted of apolar molecules (or atoms), like the noble 
gases. This leads to the conclusion that there should 
also exist intermolecular forces between pairs of ap- 
olar molecules. Whereas the aforementioned inter- 
molecular forces involving at least one polar molecule 
are classically conceivable, this third kind of van der 
Waals forces, the dispersion ones, occurring between 
two apolar, albeit polarizable, molecules, can only be 
fully understood within the framework of quantum 
mechanics. It was only in 1930 - after the devel- 
opment of quantum mechanics - that R.Eisenschitz 
and F.London 5J demonstrated for the first time how 
such a force can appear, by performing a second order 
perturbation theory on a quantum system composed 
of two atoms. Their result can be cast in the form of 
the following potential 



U(r) = 



(3) 
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where a;„io is the dominant transition frequency. An 
important result of their work is that these forces de- 



pend on the polarizabilities a of the atoms in ques- 
tion, which are related to the refractive index, and 
consequently to the electromagnetic dispersion in a 
medium composed of such atoms. The proportional- 
ity constant in this power law could then be evalu- 
ated by the two authors from fitted parameters of 
optical dispersion measurements. This evaluation, 
further developed in a paper of the same year by 
F.London alone [5], motivated London's coining such 
forces as dispersive forces in the latter article. This 
explanation successfully overturned the attempts to 
base interatomic forces between apolar molecules on 
permanent-quadrupole interaction, since it best fits 
experimental data as the van der Waals constant a 
itself. In London's words [B] (translated by the au- 
thors themselves): 

Since the van der Waals attraction, accord- 
ing to the previously accepted picture, is pro- 
portional to the square quadrupole moment, 
using the wave-mechanical model [for H2] 
one obtains, in what should he equalities, 
only 1/9 (according to Keesom), 1/67 or 
1/206 (according to Debye) of the actual 
value of the constant a of the van der Waals 
equation. 

Although the important calculation of the (— l/r^) 
power law is performed in the first article, being only 
mentioned in the second one, the latter is far more 
often cited than the former, and these forces are also 
called "London forces" . 

Dispersive forces are, then, the electromagnetic 
forces that occur between atoms or molecules pos- 
sessing no permanent electric or magnetic multipole 
whatsoever, and are due to quantum fluctuations on 
the atomic charge and current distributuons. They 
occur not only between two atoms, but also between 
macroscopic bodies, as shown for the first time in 
1932 by Lennard- Jones [7], who calculated such in- 
teraction between a polarizable atom and a perfectly 
conducting plane wall. Dispersive forces can be fur- 
ther divided into two kinds: nonretarded and re- 
tarded. London's work refers exclusively to nonre- 
tarded forces, which result when one considers light 
speed to be infinite, and the interaction instanta- 
neous. Retarded interactions, first calculated by 
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Casimir and Polder in 1948 [5], take into account 
the finiteness of interaction propagation speed, and 
in this case the dipole field of a first molecule will 
only reach a second one after a time interval of r/c, 
and the reaction field of the second molecule at the 
first one will be delayed in 2r/c. Such delay decreases 
the correlation between the fluctuating dipoles, what 
causes the retarded force to drop more rapidly with 
distance than the nonretarded one. In an atomic sys- 
tem, a characteristic time is given by the inverse of 
a dominant transition frequency ujmn , and a distance 
r is said to be long (or, equivalently, retardation ef- 
fects become relevant) when r/c > l/ujmn- It should 
be clear that nonretarded forces are a good approxi- 
mation when the molecule separation is small, which 
is the regime of validity of the London forces. 

Although Eisenschitz's and London's results were 
obtained by the use of perturbative quantum me- 
chanics, it is possible in the short-distance limit to 
estimate such forces with a much simpler method, 
known to have produced good results in calculations 
of this kind - as the atom-atom and the atom-wall 
van der Waals interactions - that goes by the name 
of fluctuating-dipoles method. This method can be 
found, for instance, in refs.[2],[l] and has also been 
shown to be useful in enabling, with few effort and re- 
quiring less background on quantum mechanics, var- 
ious discussions on dispersive forces, such as nonad- 
divity [9] or the nonretarded force between an electri- 
cally polarizable atom and a magnetically polarizable 
one [ID]. 

Our interest in this article lies on nonretarded van 
der Waals forces, and, more specifically, on the force 
between an atom and a macroscopic body. We wish 
to further develop the calculations of such forces by 
approaching a problem with curved geometry, to wit, 
the nonretarded ("London") force between an atom 
and a perfectly conducting isolated sphere. We will, 
in fact, approach this problem in two separate, inde- 
pendent ways. We first make use of the fluctuating- 
dipoles method and secondly perform this calculation 
in a way closer to the Eisenschitz and London ap- 
proach or the Lennard- Jones approach, making use of 
ordinary quantum mechanics. This naturally leads to 
a more reliable result than the first, and at the same 
time serves as correctness test for the fluctuating- 



dipoles method. 

Dispersive forces involving macroscopic bodies is of 
undeniable importance for direct experimental veri- 
fication, as seen in dJ , [12] , [13] (check also [13] and 
references therein). An aspect of van der Waals forces 
related to the interaction with macroscopic bodies is 
its nonadditivity, which leads to the fact that one 
cannot, in principle, obtain the correct van der Waals 
dispersion force in macroscopic cases by simply per- 
forming pairwise integration of the power law found 
for the atom-atom case. The reader interested in this 
feature of dispersive forces should consult [2],[5,0]. 

The modern quantum field theory explanation for 
such forces relies on the fact that there is, even 
in sourceless vacuum, a residual electromagnetic 
field whose vacuum expectation values (E(r,t)) and 
(B(r,f)) are zero, but whose fluctuations (E^), (B'^) 
do not amount to zero. This vacuum field can in- 
duce an instantaneous dipole in one polarizable atom 
(or molecule), and the field of this induced dipole, 
together with the vacuum field, induces an instanta- 
neous dipole in the second atom (molecule). It can 
thus be said that the vacuum field induces fluctuating 
dipoles in both atoms and the van der Waals disper- 
sive interaction energy corresponds to the energy of 
these two correlated zero-mean dipoles. A more de- 
tailed analysis of dispersion forces and quantization 
of the electromagnetic fleld can be found in Milonni's 
book t4J . For a pedagogical review of dispersive forces 
see, for instance, B.Holstein's paper [15] . 

As a warm-up and to establish basic concepts and 
notation, the next section is dedicated to reobtain 
the interaction between a polarizable atom and a 
conducting plane wall using the fluctuating-dipoles 
method. We then proceed in the following section 
to the interaction between an atom and a perfectly 
conducting isolated sphere, first by the fluctuating- 
dipoles method, then using ordinary quantum me- 
chanics. We end that section commenting the ob- 
tained results. In a last section we make our final 
remarks. 
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2 Calculation of the atom-wall 
London force 

It can be shown quantum-mechanically that the mul- 
tipole of the atom that contributes the most for this 
kind of interaction is the dipole (see [TB] for a demon- 
stration in the atom-atom case), thus motivating our 
picture of the atom as a dipole. Furthermore, we 
know that this dipole is not permanent, but fluctu- 
ates with a zero-mean value. The fluctuating-dipolcs 
method models the atom as constituted by a fixed 
nucleus and by an electron of charge (— e) and mass 
m. The binding force between them is taken to 
be classical-harmonic, thus leading to a harmonic- 
oscillating dipole d{t) of frequency ujq. We take, for 
simplicity, the oscillation direction to be fixed, albeit 
arbitrary, and so we have 



d(t) = -e x{t) X 



(4) 



where x{t) is position of the electron relative to the 
nucleus and x is the (fixed) unitary vector in the di- 
rection of oscillation. Since the atomic polarizability 
is defined by the expression d = aE, one can calcu- 
late the static atomic polarizability predicted by this 
model equaling the (static) force exerted by an exter- 
nal electric field E to the harmonic binding force: 



(5) 




To calculate the interaction of a dipole with a per- 
fectly conducting plane wall, we need to make use of 
the image method. The image produced when a real 
dipole d stands before a conducting plane is a dipole 
di of the same magnitude of the real one, and whose 
direction is described in figHJ 

We now proceed to calculate the equations of mo- 
tion of the electron of the polarizable atom under the 
influence of the electric field generated by the image 
dipole. Since the real dipole oscillates, so does the 
image one, and we need, then, the field created by an 
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Figure 1: Image configuration for a real dipole d in 
front of a conducting wall. The magnitudes of d and 
di are the same. Both dipoles (and the entire figure 
above) lie on the same plane. 



oscillating dipole, which is 

[3(d(r 



E(r,t) . 

3(d(r)-f)f -d(r) 



d(r)] 
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, (6) 



where r is the vector from the oscillating dipole to 
the point where the field is being evaluated and 
t* := i — is the retarded time. Fortunately, we 
are only interested in the small-distance nonretarded 
regime, which means that we can replace the retarded 
time by the time t and neglect the two last terms of 
the rhs of eq.®. The equation of motion for the 
electron of the real atom, once we project the forces 
acting on it to the direction of allowed motion (that 
is, the direction x of the dipole), becomes 



3(d, 



(2a)3 



, (7) 



where z and a are defined in fig[TJ Using the expres- 
sion of di from figUl and writing the scalar products 
as a function of the angle 0, we have 



x{t) -I- io^x{£) 



e^x{t) 1 -I- cos^ 



(8) 



But this is a simple harmonic oscillator equation 
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whose frequency is 



< 1 



6^(1 + cos^ 
47reo?Til6a'^cijQ 




We now, following the chosen method, quantize the 
system merely turning classical harmonic oscillators 
into quantum ones of same frequency. When the ex- 
ternal fields alter the frequency of the oscillator, we 
can quantum-mechanically say that its energy was 
altered too. Now, the essence of the method lies on 
identifying the zero-point energy variation as a po- 
tential energy, i.e., U := h{uj — ujq)/2. This means to 
compute the difference in energy between our system 
as it is and the corresponding existing system if there 
were no electric field and interpret that difference as 
an interaction potential. Assuming the atom to be 
isotropic, we replace cos^O by its spatial average of 
1/3. We thus find as a leading term 



U{a) 



h{uj — LUq) 
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(11) 



and this is the atom-wall dispersive potential found 
by this method. The frequency ujq artificially intro- 
duced before is identified with a dominant transition 
frequency of the atom. 

Different textbooks on basic quantum mechanics 
present the calculation of the atom-wall van der 
Waals' potential using ordinary quantum mechanics, 
as for instance, the one by Cohen- Tannoudji et al. 
|17j . One can in this fashion reobtain the Lennard- 
Jones result of 1932 [7] for the short-distance disper- 
sive interaction between a polarizable ground-state 
atom and a perfectly conducting wall, which can be 
written as 



U{a) 



ml 
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(12) 



where dx,dy,dz are the dipole component operators 
and |0) the ground state. If we assume a dominant 
transition, the result can be cast into an expression 
which equals the triple of the end result of eq. ITT]) . 



This similarity is typical for the fluctuating-dipoles 
method: it yields a result whose disagreement with 
quantum-mechanical results is only a constant fac- 
tor. All dependences on parameters as a, a, LUko are 
correctly displayed by this semiclassical method. 



Calculation of the 
sphere London force 



atom- 



We now calculate the van der Waals dispersive in- 
teraction between a perfectly conducting isolated 
sphere and a polarizable atom. We first reuse the 
fluctuating-dipoles method in this more involved ge- 
ometry, then proceed to a quantum-mechanical ap- 
proach. 



3.1 Calculation by 
dipoles method 



the fluctuating- 



AU previous considerations about the atom made in 
the last section apply, especially its polarizability and 
the interpretation of its frequency luq as a dominant 
transition frequency. We need to study the reaction 
of the sphere to the presence of a dipole, in other 
words, the images produced in this situation. This 
classical problem has a more complex geometry, but 
has been solved by F.C.Santos and A. C. Tort in [TS] . 
and the solution is described in figO We define 
Zr := R + a, where R is the radius of the sphere 
and a the minimum atom-sphere separation. There 
is one image dipole at radius Zi :— R^ / z^, a first 
image point charge qi at the same position and a sec- 
ond image point charge {—qi) at the center of the 
sphere. The magnitude of these images is signaled in 
figO Minding that, as before, forces due to the fields 
of the images must be projected onto the allowed di- 
rection, we find the following equations of motion for 
the electron of the polarizable atom 



x{t) + ljL!Qx{t) 



TTl'i'Keo 



qi'i 



1 



1 



{Zr 

3(d, 



Z^f 



(13) 
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di = 5(2(d-z)z-d) 




Figure 2: Image configuration for a real dipole d be- 
fore a conducting sphere. The values of the image 
dipole and image charges are indicated. R is the 
radius of the sphere, a, the sphere-atom separation, 

Zr := R + a and Zi :— R^ / z^ 

With the values of the images to be seen in figl21 
expressing the scalar products as function of the angle 
9 and using eq.Q to extract the factor x{t) out of 
each term, we get 

■ R cos^ ( 



xit) + U}^x{t) = 



1 



{Zr 



Z^)^ 



Replacing a ~ /[muj^) and expanding in a Taylor 
series 

-Rcos^ef 1 1 



2(47reo) 

i?3 (1 + 008^6 



(16) 



^ being a generic form to refer to terms of O [R/z!^) 
as well as terms of 0{a^^), the latter occurring only 
when a <C -R. We will assume from now on the value 
of (a^/47reo) to be smaller than unity, and so will our 
calculations have this validity constraint. 

Once again we turn classical harmonic oscillators 
into quantum ones, replace cos'^9 by 1/3 and take the 
zero-point energy difference h(uj ~ wo)/2 as a poten- 
tial. The leading term amounts to 



U 
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z^ V Zi)'^ Zj 



(17) 



and setting the expression in terms of the parameters 
R and a: 



UR{a) = - 



hujQa 
(47reo)12 
R 



AR^ 



{2R^ 



R 



{2R + afa? {R + af 



(18) 



This is again a simple harmonic oscillator equation 
whose original frequency cjo was altered to 



1 
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Rcos^e 



1 



i?3(l- 



-COS 



Z'^ {Zf Zi^"^ 



1 

1/2 



(15) 



The first term is due to the image dipole, the second 
one, to the charge qi and the third, to the charge —qi, 
located at the center of the sphere. 

3.2 Calculation via perturbative 
quantum mechanics 

(14) We now perform a more sophisticated, quantum- 
mechanical calculation of the dispersive interaction 
between a ground-state polarizable atom and a per- 
fectly conducting isolated sphere. We shall work on 
the Schrodinger picture, treat the fields classically, 
and deal with quantized atoms, seeking a leading 
term in our potential by use of perturbation theory. 

We now again depict the atom as a fluctuating 
dipole, this time in a different way. We shall state 
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that the images' fields vary too little with respect to 
position between the electron and the nucleus, so that 
we evaluate the field at only one point. This is the 
dipole approximation. However, some care must be 
taken when writing the interaction hamiltonian be- 
tween the dipole and the conducting sphere. It can 
be shown that the interaction energy of a generic con- 
figuration of the classical system formed by a dipole 
and a conducting sphere is not simply — d • E, but 
— (l/2)d • E (see the Appendix for a demonstration 
of this fact, where the energy of the configuration is 
computed as the external work to bring the dipole 
from infinite). 

Our quantized hamiltonian includes the atom 
hamiltonian Hq (kinetic term plus coulombic attrac- 
tion to the nucleus and repulsion from other elec- 
trons) plus this interaction energy: 



H = Hn 



1. 



E 



(19) 



where d is now the atom's electric dipole opera- 
tor, which equals the electron charge (— e) times the 
electron's position operator. We shall consider the 
last term as a time-independent perturbation to the 
atomic eigenf unctions. 

We need to deal once more with the images created 
in a sphere by the presence of an electric dipole. The 
classical picture of fig El still holds, and the field at 
the position of the atom can be split into two contri- 
butions, one due to the image dipole 

3(di • z)z - di 



E, 



-3(d • z)z + 6(d • z)z - 



d-2(d-z)z 



((d. 



A^eoizr - z,rz^ ' ^^^^ 

where we used the relation between d^ and d, shown 
in figUl and the other due to both image charges. 



E, 



E_ 



A'n:et){zr — Zif Attcqz^ 



1 
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{Zr Zi)^ 



(21) 



The perturbation hamiltonian W becomes 



1 



W = — d-E 
2 



((d ■ z)^-Kl2) B? {d-ifR( 1 1 



\'Kei^'l{zr — Zi)^z1 47reo2z2 



(yZir- Z-^y^ 

1 



(22) 



The first perturbative correction to the system 
energy is given by (OjW^lO) (with |0) referring to 
the atom ground state). For an isotropic atom, 
(0|d2|0) = (0|4|0) (0|d2|0). The energy correc- 
tion is then 



47reo2 
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Z^f 
1 



R 



- To ■ (23) 



We identify this energy shift as the interaction po- 
tential between an atom and a conducting sphere. If 
we also write the potential as a function of the pa- 
rameters R and a, we get 



(0|4|0) 
47reo2 



4i?3 



{2R + a)^a^ 
R 



R 



{2R + aYa? {R + af 



(24) 



And once again the three terms are due to d^, qi and 
{—Qi), respectively. 

Expressions like (0|(i^|0) can be calculated for a 
hydrogen atom, but, for the sake of comparison with 
the previous model, we shall assume the atom has a 
dominant transition frequency, say, ojmo- It has been 
shown (see [TH]) that the static polarizability of such 
an atom is 



d,|0) 



(25) 



where ujko is the transition frequency between the 
ground state |0) and any excited state |fc). The dom- 
inant transition assertion allows us to say that for 
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every k ^ m the term |(0|d2:|fc)P is negligible com- 
pared to I (0 1 I m) p. We can then replace every Uko 
by LOmO, since for k ^ m the contribution of each k to 
the sum will be negligible, and this allows us to take 
ujmo out of the summation. In a two-level atom, this 
step would be exact. In order to use the closure re- 
lation, we use that (Old^^lO) = in an isotropic atom 
to add the term |0)(0| to the sum, and wc find 



fi^„oa = 2(0|4|0} . 



(26) 



Hence, for a two-level atom, 



47reo4 



{2R + a)^a3 
R 



R 



(2i?-f a)2a2 {R- 
3.3 Discussion of the results 



(27) 



We now proceed to describe the most important fea- 
tures of the result, the potential given by eq. ([27|) (or, 
more generally, by eq.([24|)). The first striking con- 
sequence of this result is that the semiclassical one, 
eq.fUS]), differs only by a factor 3. A similar thing 
has already occurred in section 2, and also happens 
in the atom-atom (nonretarded) case: the results by 
the two methods only differ by a constant prefactor, 
and this factor does not alter the order of magnitude 
of the interaction. Our second result is clearly more 
reliable, though. 

There is a very good way to check our result with 
the known calculations on this subject. One only 
needs to take the limit i? — > oo keeping a constant, 
in which the conducting sphere would turn into a 
conducting plane wall. Wc do this directly from the 
general expression of eq. ((24|) , finding the potential 



Uooia) 



m\o) 

47reo4a3 



(28) 



which coincides with the previous known result given 
by eq. (jl2p for the dispersive interaction potential en- 
ergy between an atom and a perfectly conducting 
plane wall in the nonretarded limit. 

Another interesting limit we can take is more pe- 
culiar, in the sense that it corresponds to an inter- 
pretatively challenging physical situation: the limit 



i? ^ a, in which the sphere would turn into a so- 
called conducting point: 



4^0206 



(29) 



This is formally equivalent to the asymptotic behav- 
ior as i? —> 0, a constant. Although the physical 
interpretation of what a conducting point represents 
is rather subtle, eq. (^5)) can serve as a very useful ap- 
proximation for situations in which the conducting 
sphere is much smaller than other distances in ques- 
tion. Furthermore, there is a strong resemblance to 
the London atom-atom interaction. Since a/Aireo has 
the dimensions of a volume, it can be loosely inter- 
preted as an effective atom volume in its interaction 
with photons. London's result (eql3]) would consist 
of the product of the transition energy with both ef- 
fective volumes over r^. The corresponding volume 
of the sphere would be its real volume, and, pref- 
actors aside, eq. (f29l) also consists of the transition 
energy times the effective volumes of the atom and 
the sphere over a^. These results are, in that sense, 
equivalent to each other. 

We lastly present a graph of the obtained potential 
as a function of the distance a between the atom and 
the sphere, figlS] Using arbitrary units, we set the 

mm 



prefactor 



equal to unity, and R = 0.5, where 



47reo2 

a has the same units as R. The graph is monotonic, 
and we separate it in three terms, Ud, the contri- 
bution from the image dipole, /7+g, the contribution 
from the image charge closer to the surface, and U-q, 
the contribution from the image charge at the center 
of the sphere. One can clearly see that the term U-q 
contributes to a repulsive potential, although the to- 
tal interaction is always attractive. It is intuitive, and 
can be demonstrated in ea. ([M)) . that \U-q\ < \U+q\, 
since the charge —qi is farther from the dipole. As 
the value of i? increases, the total interaction for each 
a increases in modulus, although the images grow 
farther from the atom. This happens because with 
increasing R the graph becomes dominated by the 
term Ud, and the image dipole has a [R/ Zj-Y depen- 
dence. The terms U+q and actually decrease (in 
modulus) with increasing R. 
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Figure 3: Graph of potential between atom and 
sphere, eq.l^. We set i^MilM ^ u^ia) is in 
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arbitrary units and R and a are in the same arbi- 
trary units. We separated three contributions for the 
potential, Ud due to image dipole, U+q due to the 
closer image point charge and U-q, due to the far- 
ther image charge. 



the fluctuating-dipoles method and of the ordinary 
quantum-mechanical van der Waals calculations in 
situations with curved geometry. We performed one 
important accuracy check on our result, which is the 
limit R oo, in which the sphere turns into a plane. 
Our final result does agree with the literature in that 
limit. 

The most remarkable feature of our result is the 
calculation of the force when the sphere turns into a 
conducting point. Besides its possible approximative 
value, the ideia of a conducting point may be useful 
in more complex situations, as for instance in the 
simulation of defects in field theory. 

Our first perspective is to find the retarded disper- 
sive interaction between an atom and a conducting 
sphere. A reader might be attempted to consider the 
calculation of the long-distance, retarded interaction 
by use of the fluctuating-dipoles method or of per- 
turbative quantum mechanics without neglecting the 
radiation fields as we did after eq.®. However, the 
retardation effects require the use of a more com- 
plete description, one including the quantized fields. 
We also envision to work with dispersive interactions 
with curved geometries in general, such as rugged 
surfaces, and to generalize our results for materials 
of limited conductivity, such as dielectrics. 
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4 Final remarks 

It should be stated that all three kinds of van der 
Waals forces have many different practical and theo- 
retical applications, for instance, in condensation and 
crystallization, in structural and energetic effects in 
colloidal chemistry or biology, in the vast field of ad- 
hesion (including its applications in washing - the 
role of a detergent is to diminish the van der Waals 
forces between dirt and tissues), in its connection to 
the Casimir effect, among others. 

We have performed the calculation of the van der 
Waals nonretarded dispersive force between a po- 
larizable atom and a perfectly conducting isolated 
sphere. It constitutes an interesting application of 
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Appendix 

We now calculate the energy of a classical system 
composed of a dipole and an isolated conducting 
sphere. We shall do this computing the external work 
required to bring the dipole from infinity into its final 
position using a particularly simplifying path. Our 
result does not lack generality, though, since the end 
configuration is arbitrary and, as we know, this work 



9 



is path- independent. 

We first bring the dipole from infinity keeping it 
parallel to the y direction (see figE]), or keeping 
9 = TT /2. In this setup we only have one image, the 
dipole, of the form d, = -dR^/{R + af. The field 
generated by this image is 



The field on the real dipole is 



E = 



3(di • z)z - q.^z / I 1 



./3 



47reo V r' 



ITo - -7 (36) 



E(r') 



3(d, 



R'd[3{-y/r')r' + y] 
{R + a)347reor'3 



(i?+a)347reor'3 {R+a)^ATreor'^ 
d,_Ri / I 1 



(37) 



. ^^^^ 

where r' is the vector from the image to the point of 
evaluation of the field, and r' = \J [z — ZiY + y^- We 
shall need the field for y ^ to derivate the field in 
the next step. The force on the dipole obeys, on our 
case, 

d 

F = (d • V)E = (dy • V)E = d— E (31) 



The torque on the dipole is of the form d x E, and 
its only relevant component is on the x axis (out of 
the page). We thus need to know dyE^ — d^Ey, which 
equals 



dydz 



R^ 



R f I 1 



47reo ir'3(i?+a)3 
The work on this second step is 



(38) 



3j2 



R'd 



{R + a)347reo 



3f' , ; 



, (32) 



Wii = - J {d X F,) ■ de = - J {d X E)^de 



where (...) is a nonsingular vector quantity, and 

p _ -3d^zR^/{4TTeo) _ -3d^z R^{R + a) 
~ {zr-z^)*{R + a)^ ^ Atteo a*{2R + a)^ 



(33) 



The work done in this first step is Wi — — / F • dr 



Wi 



id^R 



2 p3 i-rf 



[R + a')da' 



« (1 + OdC 



Aneo J^d'H2R + a'Y AireoR^ ^(2 + 0^ 
-3d2 i?3 -d2^V(47reo) 



All the terms in curly brackets on eq. f38|) are taken 
out of the integral. Using that dy — d s'mO and 
dz — d cos 0, the integral to perform is 



d^ f sine' cos 6' de' = d^S^ ^ , (39) 

Jtt/2 2 2 

where in the last expression (and from now on) d^ is 
the final dipole component. 



47reo6a3(2i? + a)3 2{zr - Zi)^ {R + af 



(34) Wii 



R' 



R f I 1 



47reo2 yr'^{R+af z^ \r' 



The reader challenged by the above integral can make 
use of analytical integration softwares, such as Maple 
or Mathematica, to find its surprisingly simple result. 

We now proceed to rotate the dipole into its final 
position, that is, from 9 — 'k 12 to an arbitrary 9 in 
figlH and calculate the work done by the torque on 
the dipole. Besides the image dipole, image point 
charges qi and —qi are now present. 



Since E„ 



and 



(40) 

Zr — Zi, one can rec- 



ognize, comparing to eq. (PT|) . the last term as 
— (l/2)d • (E+q -I- E_q). Summing the first term with 
Wi from eg. pi]) and using that 



R^{d^ + dl) _ R'jMl + dl-dl) 
[R + a)3 (R + a)3 

we find that 



d-[3(d,-z)z-d,] 

(41) 



d,R 

1^ = 



di = {dzz. ~ dy-y) 



i?3 



[R+af 



(35) 



Wi + Wii = -id.[Ed,+E+g+E_,] = -^d-E (42) 
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justifying eq. p^ . Our result holds for any value of 
R, which includes the limit R — > oo, i.e., the plane 
wall. In the case of the wall, where we only have one 
image, the dipole, it is quite intuitive that there be a 
factor (1/2). One must only consider the fields' en- 
ergy density, and the fact this density is zero in half 
of the space (the wall itself) and equal as in the case 
of two real dipoles in the other half. Therefore, the 
energy of the dipole- wall configuration is, by symme- 
try, half the energy of two appropriately correlated 
real dipoles. The spheric conductor does not feature 
such symmetry, thus requiring the calculation in this 
appendix. 
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